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Motivated by tlie discovery of a plenitude of metastable vacua in a string landscape and the 
possibility of rapid tunneling between these vacua, we revisit the dynamics of a false vacuum bubble 
in a background de Sitter spacetime. We find that there exists a large parameter space that allows the 
bubble to collapse into a black hole or to form a wormhole. This may have interesting implications 
to inflationary physics. 
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I. INTRODUCTION 

There has been a long interest in the study of the dy- 
namics of vacuum bubbles in different background cos- 
mology [11-01 ■ Presumably, these bubbles may be formed 
in the course of phase transitions in the early Universe or 
arise from inhomogeneities generated by quantum fluc- 
tuations during inflation. In these studies, the bubble 
is modelled as a spherically symmetrical region of true 
or false vacuum which is separated by a domain wall 
from the exterior world such as Minkowski, de Sitter, 
or Friedmann-Robertson- Walker spacetime. The mo- 
tion of the domain wall, which is analyzed under the 
general-relativistic thin-wall approximation, character- 
izes the global structure of the spacetime. For example, 
dependent on initial conditions, an observer in the in- 
terior will describe a inflating bubble universe while an 
observer in the exterior will describe the bubble as a black 
hole. 

Recently, the possibility of a string landscape pj has 
opened a new perspective on the vacuum bubble dynam- 
ics. The string landscape is considered as a complicated 
potential energy function in a multidimensional moduli 
space with a plentitude of local minima separated by bar- 
riers of various heights. One of these minima corresponds 
to the vacuum where our Universe exists today and the 
evolution of the very early Universe is determined by the 
route that the Universe has taken in the landscape. How- 
ever, a detailed mapping of the structure of the landscape 
is still elusive. 

In the string landscape, a local minimum is a 
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metastable de Sitter vacuum with certain energy scale, 
tunneling down successively through potential energy 
barriers to vacua with lower vacuum energies via bub- 
ble nucleation. When tunneling is taking place, various 
regions of the spacetime are inflating in various de Sitter 
vacua, resulting in a very complicated global spacetime 
structure with a hierarchy of vacuum bubbles separated 
by domain walls. Typically, the height of the barrier be- 
tween the vacua is of string scale, so the tunneling rate in 
the landscape is expected to be exponentially suppressed 
and far below the expansion rate of the inflating region. 
Consequently, most part of the spacetime undergo eternal 
inflation, with some regions occasionally tunneling to a 
landscape with a flat potential that drives slow-roll infla- 
tion or to anti de Sitter vacua with negative cosmological 
constant. 

Development in the understanding of false vacuum de- 
cay in the string landscape has recently been made. It 
was found that under certain conditions the decay rates 
can be greatly enhanced even for tall barriers due to res- 
onant tunneling in the presence of multiple vacua 0-[iH ; 
due to stringy effects such as tunneling via the Dirac- 
Born-Infcld action [13 , by an effect that facilitates 
giant leaps in the landscape [ij], or by having exponen- 
tial prefactors [l3|- In light of this, it is possible that 
the tunneling rate is not smaller than the expansion rate 
of the de Sitter phase and thus a few or even copious 
vacuum bubbles may be formed within a Hubble radius. 
This may lead to bubble collisions [l^ with one of the 
domain walls near our horizon, or the formation of black 
holes and wormholes that induce new metric perturba- 
tion during slow- roll inflation [T6| . Furthermore, sam- 
pling the multiple vacua with shallow barriers by quan- 
tum fluctuations during inflation may populate domain 
walls that may affect the subsequent evolution of the Uni- 
verse [13 • AH of these, though still speculative, may leave 
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important observable signatures of the string landscape. 

In light of this, we reexamine the dynamics of a de 
Sitter bubble in a background de Sitter spacetime. The 
bubble is assumed as spherically symmetric and sepa- 
rated from the background by a thin shell. Then we use 
Israel's thin shell formalism [l^ to glue the bubble to 
the background, as was done in Refs. [lMl,i]. Although 
most of the global structure of the spacetime has been 
categorized, the parameter space that we are interested 
here has not been fully explored. The paper is organized 
as follows. In Sec. |lTl we review Israel's thin shell for- 
malism. In Sec. mil we derive the equation that governs 
the motion of the shell and show that in general the shell 
has peculiar velocity. In addition, we briefly review the 
maximally extended coordinate systems of the de Sit- 
ter spacetime and the Schwarzschild-de Sitter spacetime. 
The evolution of the bubble is discussed in Sec. IIVI and 
the numerical results are summarized in Sec. |Vl Sec- 
tion IVII is the conclusion. Here we follow the sign con- 
vention and definitions in Ref. uM. 



II. THIN SHELL FORMALISM 

To describe the behavior of the shell which divides the 
spacetime into two regions A^"*" and it is simplest 

to introduce a Gaussian normal coordinate system in the 
neighborhood of the shell. If we denote by E the time- 
like (2-|-l)-dimensional spacetime hypersurface in which 
the shell lies, we can introduce (2-t-l)-dimensional coordi- 
nates = (r, x'^ , x^) on E, where i runs from 1 to 3. Here 
T is the proper-time variable that would be measured by 
an observer comoving with the shell. In the Gaussian 
normal coordinate system, a geodesic in a neighborhood 
of E which is orthogonal to E is taken as the third spatial 
coordinate denoted by rj. Thus, the full set of coordinates 
is given by x^ = (x* , ij) . We can take x^ = 6 and x^ = cj) 
in the case of a spherically symmetric shell. Then the 
metric can be written as 

ds^ ^ ^ g^^dx^dx^ ^ (1) 

We introduce an unit vector field £,'^{x) which is normal 
to 77 = const, hypersurfaces and pointing from A^~ to 
. The induced metric on the hypersurface E can 
then be written as 

hf_iiy = Qfiiy ^flS.L'- (2) 

In the Gaussian normal coordinates, 

^'^(a;)=eM(a;) = (0,0,0,l) (3) 

The Gaussian normal coordinate system is suitable for 
the Gauss-Codazzi formalism [l^ in the neighborhood of 
E with the coordinate 77 orthogonal to the slices near E. 
The extrinsic curvature is defined by 

Kij = , (4) 



or equivalently 

K^, = = ^9,gy- . (5) 

It is evident that Kij is a symmetric three-tensor. 
The Einstein equations are 

^A"' ~ ^9tiiyR = SttCT^i. , (6) 

where is the Ricci tensor, R is the Ricci scalar, and 
T^i, is the matter energy-momentum tensor. In addition, 
the energy-momentum conservation law implies 

T^'-u = , (7) 

where the semicolon represents the four-dimensional co- 
variant derivative. Assuming that the energy-momentum 
tensor only has ^-function singularity on the shell, one 
can rewrite it in terms of Gaussian normal coordinates 
as 

Tp"- = Sf,^{x')S{r]) + (regular terms) , (8) 

where S'^jy is the surface stress-energy tensor of the shell. 
Adopting the thin shell formalism we can rewrite 
the Einstein equations and the energy-momentum con- 
servation law as 

[K/] = -8^G(5/ - ^S)S) , (9) 

. + m"] = , (10) 

and 

{K,^}Si - [T,^] ^ , (11) 

where = ^{K± + Kr^), [A] = A+ - A' for any 

quantity A, and a subscript vertical bar demotes the 
three-dimensional covariant derivative. 

By the symmetry analysis and the energy-momentum 
conservation law [2] , we have 

Sf"" = aiT)^^ - C{T){h'"' - VU") , (12) 

where = (1,0,0,0) is the four-velocity of the shell 
in the Gaussian normal coordinates. Here cr and C, are 
the surface energy density and the surface tension of the 
shell, respectively. In the case of a spherically symmetric 
shell, the induced metric on the shell can be written as 

ds^\^^^dT^ +r^{T)d9? , (13) 

where rij) is the proper circumferential radius of the 
shell. Eq. pUj) then reduces to 

a = -2(fT-C)-, (14) 
r 

where the overdot denotes a derivative with respect to 
T. Using the thin-shell approximation and the under- 
lying field-theoretical dynamics of the scalar field which 
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comprises the shell one can show that ( = a. From 
Eq. dm), it follows that & = 0. Thus, Eq. (HI]) reduces 
to 



(15) 



where ct is a constant independent of t. By combining 
Eq. jg]) with Eq. one finds 



(16) 



III. JUNCTION CONDITIONS 

For the system under consideration, we have 

_ J — (Ai/87rG')g^^ interior of the shell, , > 
^ \ -{A2/8TrG)g^,^ exterior of the sheh, 

where Ai and A2 are the cosmological constants of the 
shell interior and exterior, respectively. For convenience, 
we will consider positive values for both Ai and A2 in 
the following. At the end, we will extend the results for 
negative Ai and A2. 

^From Birkhoff theorem, we can easily write down the 
most general 0(3)-symmetric solutions to the Einstein 
equations for a region of spacetime with a non-vanishing 
cosmological constant. They are the de Sitter (dS) and 
Schwarzschild-de Sitter (SdS) spacetimes. Thus, the met- 
rics of the interior and exterior spacetime in the static 
coordinate system are 



dsl = -F±{r)df + F±-\r)dr^ + r^dil^ , (18) 



where 



F±(r) 



interior spacetime, 
1 — 2MG/r — X'2^r^ exterior spacetime. 



(19) 

Here M is an as yet undetermined parameter, xi^ = 
Ai/3, = A2/3, and dVi^ = d9^ + sin^ Odcj)"^. The cor- 
responding 6-9 components of the extrinsic curvatures of 
S, as it is embedded into the de Sitter and Schwarzschild- 
de Sitter spacetimes, can be easily calculated in the static 
coordinates by a standard recipe [l|, H, [l^ • They are 



KeeidS) = r/3, 



dS 



and 



where 



and 



KgeiSdS) = r/3sds 



/3ds-±(l-XiV+f2)i/2 



/3sds = ±(1 - 2MG/r - X2^r^ + r'^f'^ 



(20) 



(21) 



(22) 



(23) 



Substituting Eqs. and into Eq. (HH), we get the 
equation of motion of the shell. 
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r = 



MG{k^ 



X2 



where 



47rGCT > 



-1, (24) 



(25) 



and 



X4 



Xi 



X2 



(2kx2 



(26) 



It is perhaps impossible to analytically solve Eq. 
However, the asymptotic behavior of r is easily obtained. 
We have 



and 



r(T) oc exp(x-|_^r/2K) 



nT) oc I — - — 1 



00, 



0. 



(27) 



(28) 



Now we are going to show that the shell is not static 
relative to an observer using the comoving coordinates. 
Since the Schwarzschild-de Sitter spacetime approaches 
the de Sitter spacetime asymptotically, it serves our pur- 
pose by simply using its asymptotic form. For a de Sitter 
spacetime, the transformation between the static coordi- 
nates (t, r) and the comoving coordinates (f, f ) is well- 
known, given by [20| 



f = re-X*|l - ;^;2^2|-l/2 



(29) 



where we have temporarily suppressed the subscript of 
X for convenience. ^From the viewpoint of an exterior 
observer in the asymptotic region, the proper circumfer- 
ential radius r is related to the comoving radius f by 



r 

W) ' 



(30) 



where a{i) = exp(x2?) is the scale factor of the exterior 
spacetime. The proper-time rate of change of r is 



dT ) 



SdS 



r - X2r{l - X2^r-^ + f^)^/^ 
2I21 



a\l - X2^r' 



(31) 



where a is now given as a function of r. Compared with 
Eq. (|24|) . we find that in general df/dr 7^ 0, i.e., the shell 
cannot be evolving with a fixed comoving scale. Simi- 
larly, the proper-time rate of change of t is 



SdS 



(1 - X2^7'^ + r^)^/^ - X2rr 
|1-X2V2| 



(32) 



Using dr/dt — {dr / dr) / {dt / dr) and Eq. 

2kx2 



SdS 



1 



X2^ 



we have 
1/2 



(33) 
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as r — )■ cxD. On the other hand, from the viewpoint of 
an interior observer, we can repeat the above arguments 
and obtain 



dtJds 



a-' it) 



2kxi 



- Xi + X2" 



-1-1/2 



(34) 

as r — >■ Qo. Hence, the necessary condition for a sheh 
moving with a fixed comoving scale is 



\Xi 



X2 



(35) 



This means that the surface energy density of the sheU 
has to compensate the difference of the vacuum energy 
densities across the shell. Moreover, the shell cannot be 
comoving with respect to both observers simultaneously. 
Note that if k = 0, the shell will become null. 

Before solving the evolution of the bubble, we should 
address some issues regarding the topological structure 
of the spacetime. Firstly, there is a sign ambiguity in the 
expression of the extrinsic curvature (|22p and in Eq. 
as well. Moreover, in the static coordinate system, there 
exist coordinate singularities at the event horizon, which 
are indeed caused by a poor choice of coordinates. Never- 
theless, the sign of the extrinsic curvature and the phys- 
ical interpretation are manifest when it is expressed in 
terms of the Kruskal-like maximally extended coordi- 
nates. 

The Gibbons- Hawking maximally extended coordinate 
system (f 0, </>) for the de- Sitter spacetime is well 
known [l|, [21[ , and it is related to the static coordinate 
system for r < by 



1 - Xr 
1 + xr 

1 - xr ' 
1 + xr 



1/2 



cosh(xi) 



1/2 



sinh(xi) . 



The new metric is given by 

ds^ X"^(l + xr?{-dv^ + du^) -f r^d^"^ , (36) 



where r is a function of u and v given by 
1 - X?- 2 2 

= u — V . 

l + Xr 



(37) 



The metric (|36|) can be maximally extended over the en- 
tire u-v plane, subject to the constraint — < 1 . 
The coordinate system then covers the entire de Sitter 
spacetime. In terms of the Gibbons-Hawking maximally 
extended coordinate system, we can rewrite /3ds ^ 



/3ds 



-Xi ^ {I + Xirf - vii) 



(38) 



In the case of the Schwarzschild-de Sitter spacetime, 
the maximally extended coordinates was first derived by 
Bazanski and Ferrari [2^]. Unlike the Schwarzschild or 



the de Sitter spacetime, the entire Schwarzschild-de Sit- 
ter spacetime cannot be covered by a single coordinate 
patch. This arises due to the fact that there exist both 
black-hole and cosmological event horizons. 

Recall that the Schwarzschild-de Sitter metric is given 

by 



ds^ = -F{r)dt^ + F-\r)dr^ + r^dU^ 



where 



F{r) = 1 - 



2GM 



2 2 

X r 



(39) 



(40) 



The location of the horizon in the Schwarzschild-de Sitter 
spacetime can be found by solving the equation 



(41) 



F{r) = . 



We define an extremal mass by Moxt = l/(v27xG'), 
which turns out to be a relevant mass scale in our analy- 
sis. When M > Moxt, Eq- (HB) has no positive root and 
hence there does not exist any horizon. In the case of 
< M < Mcxt, there are two positive roots, th and rc 
(rjj < rc), which satisfy 



and 



rn^ + rare + rc^ = X ^ 



r-arcirii + rc) = 2Afx" 



(42) 



(43) 



When M — A/oxt, two positive roots become identical. 
We define the quantities and dc for future use by 



1 - 3x2r^ 



5c ^ — 



rc 



1 



(44) 



The coordinate patch which can be maximally extended 
to cover the entire neighborhood of the black- hole horizon 
is given by 



(r - rH)(r + th + rc)" 



(rc - ry 
(r - rH)(r + rn + rc)^"^ 



(rc - ry 



cosh 



sinh 



t 



2fe 
t 



2Sh 



where rjj < r < rc and x = (5c/<5h- The new metric is 
given by 



ds^ = Ax^5l{rc~r)^+''{r + ru + rcf~''i-dv^ + du^) 
+r^dn^ , (45) 

where r is a function of u and v given by 

(r - rH)(r + rn + rc)^"^ _ 2 2 



(rc - ry 



(46) 



The metric (j45p can be maximally extended over the en- 
tire u-v plane, subject to the constraint 



v"^ > -rnr^'-^' {r-n + rc) 



x-l 



(47) 
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The coordinate system then covers the entire neighbor- 
hood of the black-hole horizon in the Schwarzschild-de 
Sitter spacetime. 

Similarly, the coordinate patch which can be maxi- 
mally extended to cover the entire neighborhood of the 
cosmological horizon is given by 



[re - r) 



(rc - r) 



cosh 



sinh 



t 



for rn < r < rc. The new metric is given by 



2Sc 
t 



25c 



+r^dn^ , (48) 



where r is a function of u and v given by 



rc 



(r - rn)- (r + r^ + rc) 



2 2 

= U — V 



(49) 



The metric (|48p can be maximally extended over the en- 
tire u-v plane, subject to the constraint 



> -1 



(50) 



The coordinate system then covers the entire neighbor- 
hood of the cosmological horizon in the Schwarzschild-de 
Sitter spacetime. 

In terms of the Bazanski-Fcrrari maximally extended 
coordinate system discussed above, we can rewrite /3sdS 
as 



/3sds 



2X2^&(''C - r)^+^(r -I- th -I- rc)^ ""{uv - vii) 
for black hole horizon, 

—2x2'^5c{r — m)^^^ {r + r-a + rcf~i{uv — vii) 
for cosmological horizon. 

(51) 



IV. EVOLUTION OF THE BUBBLE 

To discuss the evolution of the bubble and the cor- 
responding spacetime structure, we now introduce the 
dimensionless variables 



.3 ' 



and 



where 



2k ' 



2MG 



'0 = T" 



(52) 



(53) 



(54) 



Eq. ([24)1 can be rewritten as [1| 

2 



where 



E 



4k^ 



(2GM)2/3x+8/3 ' 



and 



^^^2 + 2{k^ + X2'-Xi^)/x+^ 



(55) 



(56) 



(57) 



(58) 



Note that < 7^ < 4. The equation of motion is equiva- 
lent to that of a classical particle moving in one dimension 
subject to the potential V{z). The asymptotic behaviors 
of the potential V{z) are V{z) ~ —1/z'^ for small z and 



V{z) 



for large z. One can show that d^V/dz"^ < 



for all z and V{z) has one maximum at Zm, where 

^ 1{[8 + (1 - ^^2^/^ - (1 - 7V2)} . (59) 



The maximum value of the potential is given by 



(60) 



In terms of the new variables, we can also rewrite /3ds 
and /3sdS as 



dS 



(^) 



1 



PSdS — —5— 

\ r^K J z^ 



^^+Xl^-X2' 

x+^ 



(61) 
(62) 



If < xi^ — X2^, then /3ds will change sign at z = Zds, 
where 



^dS 



X^ 



Xl^ - X2^ - 



1/3 



(63) 



If > X2^ — Xi^i then /3sds will change sign at z = zsds, 
where 



^SdS 



X+ 



+ Xl^ - X2^ 



1/3 



In case of both /?ds and /3sds changing sign, since 

1/3 



^SdS 



= 1 



2k^ 



Xi^ - X2^ 



> 1 , 



(64) 



(65) 



it follows that Zds > zsds and the equality holds for k 
0. 
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The location of the horizon in the Schwarzschild-de 
Sitter spacetime is determined by the equation, 



2GM , o 
1 X2'r2 ^ 0^ 



(66) 



and in terms of the dimensionless variables, it can be 
written as 



= V{z) + 



1 



,(67) 



This curve is tangent to V{z) at z = zsds- Similarly, the 
location of the horizon in the inner de Sitter spacetime 
is given by the equation, 1 — Xi^''^ — 0, which can be 
written as 

E ^ - (1^^ 



z 



(68) 



This curve is tangent to V{z) at z = Zds- 



V. NUMERICAL RESULTS 

The system has four mass scales: Xii X2, and M, 
which are treated as free parameters in the present con- 
sideration. In the context of string landscape, xi and 
X2 are given by the cosmological constants of the local 
minima, and k, can be calculated as long as the profile 
of the potential barrier between the minima is known. 
The initial conditions of the nucleated bubble determine 
the value of M . So far, we have assumed positive cosmo- 
logical constants. For negative ones, the derivations in 
the previous section still hold after making the change: 
X^ ~X^- A general discussion of the system has been 
given in Ref. but there lacks details about the pa- 
rameter space that we are concerned here. In Ref. 61, 
the authors considered Xi > X2 > only. In Ref. |Jj, 
only a small k was considered in the context of quantum 
fluctuations during inflation. In the following, we will 
follow the methodology in Ref. f3| that is suited to our 
purposes. 

The plots of the curves of V{z) in Eq. ((56|) and the 
event horizons in Eqs. (|67p and (|68|) are shown in Fig. [T] 
and Fig. [21 for a representative case using xi = 3.0, X2 = 
2.9, and k, = 0.7375, in unit of an arbitrary mass scale. 
Also drawn are four horizontal lines each of which denotes 
the value of E corresponding to the Schwarzschild mass 
M given in Eq. ([57]) . Mext is the extremal mass. A/max — 
V{zm), Msds = ^(^Sds), and Mds = V{zds)- 

Figure [3] shows the global geometry of the spacetime 
represented by Penrose diagrams. In each category, the 
left side is the de Sitter spacetime and the right side is 



the Schwarzschild-de Sitter spacetime. The dashed or 
solid curves with arrows denote the trajectories of the 
bubble wall as seen by the observers on both sides. The 
Schwarzschild-de Sitter spacetime in the first three cat- 
egories on the left column does not have event horizons. 
For examples, the dashed line in the diagram labelled by 
R6, R7a, R7b, R9a, R9b, or i?10 denotes the formation 
of a black hole. The dashed line in the diagram labelled 
by i?5a, R5b, R8a, or R8b denotes the formation of a 
wormhole. 

We then identify the domain of the dimensionless pa- 
rameter space (2MGx2, X1/X2) that fall into the same 
category, for three representative values of k given 
by = Ixi^ - X2^|, = 10*1X1^ - X2^|, and 
— lO^^lxi^ — X2^|. Given a set of values of 
(2MG'x2, X1/X2), both V{z) and E are fixed and the 
motion of the bubble wall is determined. The phase- 
space diagrams are drawn in Fig. 01 Fig. [SI and Fig. [51 
respectively, in which the boundaries separating different 
phases correspond to the masses defined in Fig. [H 



VI. CONCLUSION 

The global structure of the spacetime describing a vac- 
uum bubble in a background vacuum has been studied 
using the general-relativistic thin- wall formalism. In this 
work, the dynamics of the bubble wall is controlled by 
the cosmological constants on both sides and the surface 
energy density of the wall. There exists a large parameter 
space that allows the bubble to collapse into a black hole 
or to form a wormhole. The condition for the existence 
of this parameter space should be easily satisfied in the 
context of string landscape. To investigate the formation 
rate of black holes and wormholes in the string landscape 
is an interesting problem. We may speculate that while 
slow-roll inflation is taking place along a flat potential 
in the string landscape, black holes or wormholes are 
being produced by bubble nucleation in quantum tun- 
neling from the flat plateau to neighboring or even far 
away local minima with lower vacuum energies through 
the landscape barriers, or perhaps by quantum fluctuat- 
ing to nearby local minima through shallow barriers. If 
the formation rate is significant compared to the expan- 
sion rate of the inflation, this may have relevant effects 
on inflationary physics and may leave an imprint in the 
metric perturbation. Overall, this work is a theoretical 
study of the dynamics of vacuum bubbles. It itself is an 
interesting problem in general relativity, while being an 
effort towards exploring the observable signatures of the 
string landscape scenario. 
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FIG. 1: Graph of V{z) (sohd line), de Sitter (dashed line) and 
Schwarzschild-de Sitter (dotted line) horizon lines for xi ~ 
3.0, X2 ~ 2.9, and k = 0.7375 (in unit of a mass scale). Long- 
dashed horizontal lines denote the values of E corresponding 
to Mcxt, Afmax, A^sds, and Mds, respectively. 
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FIG. 2: Blow-up graph of the upper-left region around the 
peak of the potential in Fig. 1. 
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FIG. 3: Penrose diagrams representing the motion of the bub- 
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classified in the same category with an additional suffix. 
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